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The SuperKamiokande atmospheric neutrino measurements leave substantial room for nonstan- 
dard interactions (NSI) of neutrinos with matter in the — Vt sector. Large values of the NSI 
couplings are accommodated if the vacuum oscillation parameters are changed from their standard 
values. Short and medium baseline neutrino beams can break this degeneracy by measuring the 
true vacuum oscillation parameters with the disappearance mode, for which the matter effects 
are negligible or subdominant. These experiments can also search for the v^. — Vt flavor changing 
effects directly, by looking for — conversion caused by the intervening matter. We discuss 
both of these methods for the case of MINOS. We find that, while the present MINOS data on 
disappearance induce only minor changes on the constraints on the NSI parameters, the situation 
will improve markedly with the planned increase of the statistics by an order of magnitude. In 
that case, the precision will be enough to distinguish certain presently allowed NSI scenarios from 
the no-NSI case. NSI per quark of about 10% the size of the standard weak interaction could give 
Su Vfj. — conversion probability of the order ~ 10~^, measurable by MINOS in the same high 
statistics scenario. In this v^^ — channel, the small effects of NSI could be comparable or larger 
than the vacuum contribution of the small angle The expected 613 bound at MINOS should be 
more properly interpreted as a bound in the ^13 - NSI parameter space. 

PACS numbers; 13.15.+g, 14.60.Lm, 14.60.Pq 



I. INTRODUCTION 

Neutrinos offer certain unique advantages for the 
search of new forces beyond the Standard Model. These 
advantages stem from the fact that neutrinos are sensi- 
tive only to the weakest known forces, the weak interac- 
tion and gravity, and therefore offer a "cleaner" ground - 
with respect to other particles - to look for the subdom- 
inant effects of the new, nonstandard interactions (NSI) . 
The abundance of neutrinos from natural sources (chiefly 
the Sun and the interaction of cosmic rays in the Earth's 
atmosphere) adds further motivation. 

The very same reason that makes them an interest- 
ing probe, the weakness of their interactions, also made 
the neutrinos historically difficult to study experimen- 
tally. The result is that to this day some of the neutrino 
interactions are very poorly known, with uncertainties 
comparable to, or even exceeding, the standard model 
predictions. In the case of neutrino-matter interactions 
the measurements of neutrino cross sections at CHARM 
and NuTeV P, Q constrain chiefly the coupling of the 
muon neutrino, while leaving the I'e — Vt channels largely 
unconstrained (see for a comprehensive discussion). 

Measurements of neutrino matter effects at neutrino 
oscillation experiments offer increased sensitivity to neu- 
tral current NSI in the e — r sector, down to tens of 
per cent of the weak interaction strength [E IE 13 ■ Still, 
however, certain directions in the parameter space exist 
where large NSI are allowed, because their effect on the 
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oscillations of solar and atmospheric neutrinos is degener- 
ate with that of the mixings and mass square differences 

Logically, the next step towards stronger constraints 
- or discovery - of neutrino NSI is to break these de- 
generacies. One way to do so is to perform precision 
measurements of the neutrino mixings and mass squared 
differences in conditions where matter effects are negligi- 
ble or subdominant. Such conditions of "quasi- vacuum" 
are realized in the v,, — oscillation channel of the 
beams used by K2K and MINOS [Ulil. There, 
the large — (vacuum) mixing, ^23, induces oscil- 
lations effects of order unity. These dominate over the 
matter effects, which are at most of order ^ 10~^, due 
to the small matter width along the beam trajectory. 

In a recent paper |2 we showed how the data from 
K2K, in combination with the SuperKamiokande atmo- 
spheric neutrino measurements, restrict the allowed re- 
gion of the NSI couplings. The conclusion was that, while 
improving with respect to the atmospheric data only, the 
combined data still allow couplings of order unity (rela- 
tive to the Fermi constant, Gp). 

Here we discuss how the situation changes in light of 
the recently reported search for v^^ disappearance at MI- 
NOS 14]. We also illustrate the increased potential that 
MINOS will have with higher statistics. We get that with 
the present MINOS data the presence of NSI is neither 
revealed nor excluded. The allowed region of the NSI pa- 
rameters is changed only slightly. Qualitative improve- 
ments are expected, however, with the planned increase 
of statistics by an order of magnitude. The precision in 
that case would be such that in certain NSI scenarios 
(with the NSI couplings on both u and d quarks of the 
order of 10% of the weak interaction), the point of zero 
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NSI would be allowed only with confidence level higher 
than 99%. 

A second way to search for the flavor-ch anging Ve — Vj 
NSI is to look for the — conversion [T^ [T^, [13) 
[Tgj . In this channel, the vacuum contribution to the os- 
cillations of beam neutrinos is due to the small mixing 
angle 6*13. Considering that sin^ 26'i3 ^ 0.14 P0II2H we 
see that at K2K and MINOS the NSI effects could be 
of the same order or larger than the vacuum ones, and 
therefore they could be visible in spite of their smallness 
on absolute scale. In this paper we illustrate this point 
for the specific case of MINOS, and show that a conver- 
sion probability of ~ 0.05, within the reach of MINOS, 
could be due to NSI within the currently allowed region, 
with couplings per quark of the order of 10% of the weak 
interaction. 

Finally, although not discussed here, another impor- 
tant method of searching for new tau neutrino inter- 
actions at MINOS is by measuring the neutral current 
event rate in the far detector, where a large fraction of 
events should be due to tau neutrinos generated by oscil- 
lations. The combinations of the NSI parameters probed 
in this method do not uniquely map to those probed by 
the matter effects studied here lOj, hence providing a 
complementary probe of NSI. 

In Sec^Jwe describe the parametrization of the non- 
standard interactions being probed, list various con- 
straints on these interactions and give the resulting os- 
cillation Hamiltonian. We also review the main results 
of our earlier study of NSI effects on atmospheric neutri- 
nos (in Sect. Ill B|) . In Sect. IIIll we study the combined 
NSI sensitivity of the Vfj^-^Vr mode at MINOS and atmo- 
spheric and K2K data. Sec. II VI is devoted to the v^j^^Ve 
appearance channel at MINOS. Finally, Sec. summa- 
rizes our findings. 



II. NEUTRINO OSCILLATIONS AND NSI IN 
THE i/e vr SECTOR 

A. Effective Lagrangian, bounds, oscillation 
Hamiltonian 

The neutrino oscillation Hamiltonian in vacuum is 
dominated by the atmospheric mass splitting. In the 
fiavor basis i^e , , , this Hamiltonian is given by 
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where A = Am^/(4£'^), is the neutrino energy and 
Am^ is the largest mass square difference in the neutrino 
spectrum (often called Ato^j)- Corrections to H^ac due 
to the smaller, "solar" mass splitting and the remaining 
two mixing angles produce only subdominant effects on 
the v^-^Vt oscillations (see sec. IIII|) and so will be omit- 
ted for the moment. See, e.g., jlOl] for the full Hamilto- 
nian including all these terms. 



In matter, the neutrino propagation is affected by re- 
fraction (coherent forward scattering). We parameterize 
the corresponding part of the Hamiltonian as: 
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(up to an irrelevant constant). Here rie is the number 
density of electrons in the medium and the epsilon terms 
represent the contribution of the NSI. 

The physical origin of the NSI terms in Eq. ^ can 
be the exchange of a new heavy vector or scalar par- 
ticle, described by the effective low-energy four-fermion 
Lagrangian 

L^S' = -2V2Gp{9^jpiyp)ieif,''fLYh+eil''fRYfB) 
+ h.c. (3) 

Here e{/p^ {^{//' ) denotes the strength of the NSI between 
the neutrinos v of flavors a and /3 and the left-handed 
(right-handed) components of the fermions / and /. The 
epsilons in Eq. Q are the sum of the contributions from 
electrons (e*^), up quarks (e"), and down quarks (e"*) in 
matter: = J2f=u,d.e^ip'^f /"^e. They are thus nor- 
malized per electron. In turn, e^^ = e^^ + e^^ and 

= ^afi ■ Notice that the matter effects are sensitive 
only to the interactions that preserve the flavor of the 
background fermion / (required by coherence |22| |) and, 
furthermore, only to the vector part of that interaction. 

Neutrino scattering tests, like those of NuTcV 2] and 
CHARM Jj, mainly constrain the NSI couplings of the 



muon neutrino, e.g. 
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10 ^. The limits they place on Cee, eer, and Crr are rather 
loose, e. g., le™-^! < 3, -0.4 < e™^ < 0.7, |e™'| < 
0.5, |e^g| < 0.5 3. Stronger constraints exist on the 
corresponding interactions involving the charged leptons. 
Those, however, are model-dependent and do not apply 
if the NSI come from the underlying operators containing 
the Higgs fields Q. Here we are interested only in direct 
bounds that follow from experiments. 

Given the above bounds we will set e^^ and e^^ to zero 
in our analysis. Furthermore, for simplicity we will also 
set e^r to zero. The earlier analyses of the atmospheric 
neutrino data have indicated that this parameter is 
quite constrained (e^T< lO^^^ — 10^^). Here, as in our 
previous works, we have a three-dimensional NSI param- 
eter space, spanned by e^e, Cer, and err- 



B. Atmospheric neutrinos and NSI 

The atmospheric neutrino flux is comprised of both 
neutrinos and antineutrinos, initially produced in the 
electron and muon flavor states, with the energy spec- 
trum spanning roughly five orders of magnitude, from 
10-1 to 10"* GeV lij. The data can be divided into a 
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high energy sample, E^, > 10 GeV, where only the muon 
neutrino flux is measured (through-going and stopping 
muons sample), and a lower energy one, where both the 
muon and electron neutrino data are available. For the 
first, the data show disappearance of with large sup- 
pression at lower energies/longer baselines and smaller 
suppression in the opposite limits; the second sample is 
consistent with no oscillations of and disappearance 
with large amplitude. This means that the oscillate 
into Vr at low energy, and that the data are well de- 
scribed by vacuum oscillations with large mixing at all 
energies. 

One way to fit these data is to require that the NSI cou- 
plings of the tau neutrino are negligible at all energies. 
For example, if Cer^ 0, we would need err Gpne S5 
AmV(2-Bo) at ^ 20 - 30 GeV (the highest energy 
where an oscillation maximum occurs for neutrinos trav- 
eling along the diameter of the Earth), to ensure that 
oscillates into Vr with unsuppressed amplitude in the 
high energy data sample. Substituting numerical values, 
this translates into the bound | err I ^ 0.2 [||. 

In we pointed out another possibility: one of the 
two (non-zero) eigenvalues of Hmat is small with respect 
to the vacuum terms, so that at high energy can oscil- 
late into the corresponding eigenstate, Vr' , without sup- 
pression. Specializing to the case which is continuously 
connected to the standard — v-t oscillation scenario 
(i.e., Vt' becomes Vr in the limit of no NSI), this gives: 



\\+ eee + err - V (1+ ^ee ^ Crr )2 + 4| e^r P| < 0.4 . 

(4) 

Although at high energy f^'s oscillate into a mixture of 
Vt and i^e, and not purely into i/r, this is not observable, 
given the absence of Vf, data in the high energy sample. 
At lower energy, where the kinetic terms of the Hamilto- 
nian dominate, the Vf^ flux is unoscillated, as the muon 
neutrinos oscillate dominantly into Vt- 

The condition Q describes the allowed region of the 
epsilons; it has the shape of a parabola centered on the 
curve 
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It is also independent of the phase/sign of t^T- 

The conversion of z^^ at high energy mimics vacuum 
oscillations with an effective mass splitting, Am^, and 
an effective mixing angle, 6'm, given by |9lll0l|: 

Am^ = Am2 [(c2e(l + c^) - 4)^/4 + {s^eCfif\ , 
tan 26l„j = 2s2eCf3/ (c2e(l + c^) - s|) . (6) 

They depend on the NSI only through the angle /3, de- 
fined as 



tan2/3 = 2| |/(1+ 



Err) ■ 
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If NSI are present but ignored in the data analysis, the 
highest energy neutrino data would give 0m, and Am^ 
instead of the vacuum parameters. Those can be inferred 



from Eq. ©. Taking the experimentally favored value 
0m = ''■/4, we find: 
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We see that the vacuum mass splitting is larger than the 
measured one, Am^ , and the vacuum mixing 9 is smaller 
than maximal {9 < 7r/4). 

The vacuum parameters are constrained by data from 
lower energy atmospheric neutrinos, E ^ I GeV, for 
which the vacuum term of the Hamiltonian dominates. 



Calling 9min and Af 



the smallest mixing and the 



largest mass splitting allowed by those, we have two con- 
straints on (3 [lU^ : 



COS^ (3 ^ tan^ 9min , cos^ 



2At 



Am? 



,(9) 



that follow from Eq. These conditions truncate the 
allowed region Q , as shown in J^] . 

Let us emphasize that the analytics in Eqs. Q-® 
is for guidance only, since it relies on the over-restrictive 
condition 101 and on the approximation of negligible vac- 
uum terms at i? ^ 10 GeV. The corrections due to the 
vacuum terms to the conditions Q and @ introduce a 
dependence on the mass hierarchy, i.e. on the sign of 
Am^, with the inverted hierarchy (Ato^< 0) giving a 
larger allowed region of the NSI. We refer to ref. jT9| for 
details on this and on the corrections due to the mixing 
angle ^13, which break the symmetry in the phase of ter- 



III. //^ - i^r OSCILLATIONS AT MINOS 
A. General considerations for K2K and MINOS 

For the setups of K2K and MINOS the neutrino os- 
cillations in the — 1^^ channel are dominated by the 
vacuum contribution, with the matter effects playing only 
a subdominant role. This is proven by the fact that both 
beams have baselines much shorter than the characteris- 
tic length (or equivalently, the width, see which de- 
scribes refraction effects, Iref = {V2GFne)-^ ~ 1.9-103 
Km for the average density of the continental crust, 
p = 2.7 g • cm~3. Specifically, the K2K baseline is 250 
Km long, only about a tenth of the refraction length, 
and this ensures that matter effects at K2K are negli- 
gible. MINOS, with its 735 Km long baseline, covers 
about 40% of Iref, and therefore is expected to see mat- 
ter effects at a subdominant, while not entirely negligible, 
level. For simplicity here we consider the zeroth order de- 
scription of both K2K and MINOS in terms of vacuum 
oscillations only (and comment on the small corrections 
due to refraction effects at MINOS later, see sec. IIII C|l . 

In this approximation MINOS and K2K measure the 
oscillation parameters in vacuum, 9 and Am^. The differ- 
ence between 9, Am^ and the corresponding parameters 
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inferred from the atmospheric neutrino data tells us how 
much room there is for NSI, through Eqs. ((SJl-lini- 

The present MINOS data (a similar description holds 
for K2K) can be fit by a range of Am^ and 9, such 
that larger values of Am^ come with smaller values of 
9. Compared to the corresponding region from the Su- 
perKamiokande atmospheric analysis, the MINOS re- 
gion is relatively narrow in Am^, Am^ = (2.5 — 
3.65) -10-3 eV2 (at 68% confidence level, C.L.) but 
extended in 6*, with sin^ 29 ~ 0.73 - 1 (68% C.L.). The 
MINOS constraint on Am^ also improves slightly on the 
analogous restriction given by K2K. 

We notice that the region currently allowed by MINOS 
is extended in the direction (smaller 9 comes with larger 
Am^) obtained from the NSI analysis of the atmospheric 
data. Due to this accidental degeneracy, we expect that 
the present-day MINOS results leave substantial room for 
the NSI, with only a loose upper bound on (3 according 
to Eq. ©. 



B. The data analysis: MINOS, SuperKamiokande, 

K2K 

We have carried out a combined analysis of the at- 
mospheric data from SuperKamiokande and i^^ Vt 
oscillation data from K2K and MINOS. Our study has 
two parts. First, we have investigated the impact of the 
recently released MINOS data on the existing NSI con- 
straints. Second, to illustrate the expected future sen- 
sitivity of MINOS, we have simulated a future MINOS 
signal with projected higher statistics. For this simula- 
tion, we considered two scenarios: one which would give 
a signal of NSI, and another that would instead give a 
constraint on them. 

The starting point in both cases are the results from 
our NSI analysis of the atmospheric neutrino data. This 
analysis was done in and 10| and involved a fit to the 
complete 1489-day charged current Super-Kamiokande 
phase I data set |2^ |2^ , including the e-like and /i-like 
data samples of sub- and multi-GeV contained events as 
well as the stopping and through-going upgoing muon 
data events. The details of the code, provided by 
Michele Maltoni, are given in U3 have 

scanned a five-dimensional parameter space (eee, Cer, frr 
, Am23, sin^ ^23) and obtained an allowed region in this 
space. The output of that scan - a set of values on a 
five-dimensional grid - was used here. 

The output was combined with a array from the 
K2K vacuum oscillation analysis. The latter was pro- 
vided by the K2K collaboration in tabular form, and 
refers to the published analysis of Ref. T?] (Fig. 4 there). 
Since that analysis is practically independent of the mat- 
ter interactions, the same array was added to each 
(Arnig, sin^ ^23) plane. For details, see again [Tol |. 

Finally, we have performed a similar fit to the released 
MINOS data. The data used were the number of ob- 
served events in each MINOS energy bin versus the cor- 



responding predictions without oscillations as given in 
|l4j . Our analysis follows the one in as closely as 
possible. The resulting allowed region and the best-fit 
point agree reasonably well with those in the MINOS 
analysis (iJ^J, slide 60). The output of this fit was then 
combined with those of the SuperKamiokande and K2K 
analyses, giving a five-dimensional allowed region in the 
space (e Am^g, sin^ 6*23). The combined results 

were then marginalized, either over the NSI or over the 
oscillation parameters. 

For the second part, we have simulated the data with 
increased statistics that might be expected in the future 
at MINOS. We have used statistics corresponding to 25 x 
10^° protons-on-target (p.o.t.), as quoted by the MINOS 
collaboration (see, e.g., ^J), henceforth referred to as 
"MINOS-high" sample. We have also checked how our 
results change with changing the statistics to 16 x 10^° 
p.o.t. 

We have simulated future data by taking the histogram 
of the expected number of events in each energy bin pub- 
lished by the MINOS collaboration, scaling it up to the 
appropriate statistics, applying oscillations to it accord- 
ing to the set of parameters chosen and finally simulating 
statistical fluctuations in the data using Poisson statis- 
tics. It is hoped that this simplified procedure still cap- 
tures the relevant effects. The most accurate way would 
of course be to perform a full Monte-Carlo simulation of 
the detector, which could be done by the MINOS collab- 
oration. 

The simulated data were then used as an input to 
the combined SuperKamiokande-K2K-MINOS analysis 
described before. 



C. Results with present MINOS data 

The results of our analysis of the existing data are pre- 
sented in Figs. nOand Table U 

Fig. n shows the effect of adding the MINOS results 
on the allowed region in the space of sin^ 9 — Am^ . The 
contours show the case of the standard interactions, while 
the filled regions were found after marginalizing over the 
NSI (see the caption for details). The Figure also shows 
the direction in Eq. ||SJ), along which svy? 9 and Ato^ 
would be found for fixed sin^ 9^, and Am^ [t^m^ = 
2.1 • 10~^ eV^ and 9^. = 7r/4 inotivated by the best fit 
to the atmospheric data only (Iq ) and NSI varying along 
the parabola 10. 

For the case of no NSI we see that both with and 
without MINOS the best fit point lies at maximal mix- 
ing, with the interval sin^ 9 ~ 0.3 — 0.7 allowed at 
3(7 confidence level. These bounds on the mixing are 
due to the atmospheric data; the main effect of adding 
the MINOS data is to further restrict Am^: Am^ = 
(1.8 - 3.4) -10-3 eV2 at 3cr C.L. 

Both with and without MINOS, in the case with NSI 
the allowed region extends along the direction in Eq. 
(ISJ, allowing smaller mixing and larger Am^ with re- 
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spect to the standard case. Besides the restriction on 
Am^ ah'eady observed, adding the MINOS data has 
the effect to move the best fit point away from maxi- 
mal mixing. The new best fit point has sin^ 6 = 0.36 
and Am^= 2.88 -lO^"^ eV^, and it hes along the curve 
as it is plotted in the figure. This happens be- 
cause the MINOS-only best fit point - sin^ 9 = 0.33 and 
Am^= 3.05 -10"^ eV^ - hes near this curve and the 
atmospheric is very flat along this direction, once the 
NSI are included. It has to be stressed, however, that 
there is only an insignificant difference between the best 
fit point and, say, the minimum of along the direc- 
tion 9 = 7r/4. The latter point is within the 68% C.L. 
contour. 

The shift of the best fit point in the space of oscilla- 
tion parameters indicates that with MINOS non-zero NSI 
are favored, even though in a statistically insignificant 
way. To check this, we have marginalized the likelihood 
over the oscillation parameters (with fixed mass hierar- 
chy) to obtain the in the three dimensional space of 
the epsilons, x^(£ee, Ett, Cer)- We find that for the in- 
verted mass hierarchy the minimum of this marginalized 
X^ is at the point (cee, Crr, Eer) = (0.6, 0.933, 0.615). This 
point is practically degenerate with others that satisfy 
the condition and correspond to vacuum oscillation 
parameters (Eq. 0) close to the MINOS-only best fit 
point. Again, we find that the difference in x^ between 
the best fit point and the point with zero NSI is minimal 
and has no statistical significance. 

In Fig. [T] the effects of matter on the MINOS sig- 
nal were not taken into account for simplicity. We have 
checked that their inclusion leaves the results unchanged 
qualitatively, with only a change of the order of 10% in 
the extent of the upper left part of the allowed region 
(for the NSI case, shaded areas). 

Fig. 13 shows sections of the 3D allowed region along 
the plane eee= for normal and inverted mass hierarchy. 
The contours correspond to intervals of x^ with respect 
to the absolute minimum given above and represent the 
regions allowed at 95%, 99% and 3a confidence levels. 
They show only minor changes with respect to the same 
contours without MINOS [3- For illustration in Fig. H 
we show the points of minimum of x^ along the plane 
eee= 0. Those appear to be away from the origin, con- 
firming what has been discussed about the shift of the 
best fit to non-zero NSI. 

Table m gives further details on our results. It presents 
information on sections of the 3D allowed region along 
several planes of fixed e^e, thus generalizing what is 
shown in Fig. [21 



D. Expected future sensitivity of MINOS 

How will the constraints on NSI improve with the 
increase of the statistics at MINOS? We consider the 
MINOS-high setup for two examples of NSI and "true" 
oscillation parameters: (i) no NSI, sin^ 9 = 0.5 and 
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FIG. 1: Regions in the space Am^ —sva^O allowed by the 
global fit before (upper panel) and after (lower panel) the 
MINOS results, with purely standard interactions (contours) 
and with NSI (filled areas). For both cases we plot the re- 
gions allowed at 95%, 99% and Scr confidence levels for 2 de- 
grees of freedom. We have marginalized also over the sign 
of Am^ and took —1 <eee< 1.6, motivated by one of the 
accelerator bounds (see~|3). The dashed line represents the 
possible positions of the best fit point of the atmospheric data 
only in the space of the vacuum parameters for fixed 9m and 
Am^, and varying NSI along the parabola 0. We used 
Am^ = 2.1 • 10~^ eV^ and Om = 7 r/4 , motivated by the best 
fit of atmospheric neutrinos alone |2a |. 
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-1 


<0.3 


[-0.4,0.6] 
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-0.5 


<0.75 


[-0.15,1.2] 


0.15, 0.1 


0. 


<1.65 


[-0.15,2.6] 


0.438, 0.255 


0.3 


<2.1 


[-0.15,3.4] 


0.67, 0.413 


0.6 


<2.5 


[-0.15,3.9] 


0.933, 0.615 


0.9 


<2.9 


[-0.15,4.5] 


1.019, 0.615 



TABLE I: Some relevant numbers describing the 3D allowed 
region of the NSI couplings for inverted mass hierarchy. We 
give the intervals of e^T and of err enclosed by the 99% C.L. 
region in planes of fixed Cee, and the values of eeT,^TT that 
minimize the in the same planes. 
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FIG. 2: Two sections of the allowed region in the three di- 
mensional space of the NSI couplings along the plane eee= 0. 
The contours refer to 95%, 99% and 3a" confidence levels. The 
points of minimum of the 'x in this plane are also shown. The 
upper (lower) panel refers to normal (inverted) mass hierar- 
chy. The eee< regions are symmetric with respect to those 
shown. 



,2 / 



2.7-10-3eV2, and (ii) eee= 0, e^^= 0.81, eer= 0.9, 
sm-y = 0.27 and Am^ = 3.1 • lO'^eV^. Fig. El gives 
the results of our fit to simulated MINOS data. In both 
cases (i) and (ii), the fits were done under the assump- 
tion of no NSI [33 . The aim was to see if a large NSI 
scenario would result in a strong tension between the MI- 
NOS data and the atmospheric neutrino fits, if both were 
analyzed assuming standard matter interactions. Indeed, 
this happens for scenario (ii): the high statistics MINOS 
allowed region has no overlap with the current region 
from MINOS+atmospheric-|-K2K. This way, the MINOS 
measurement would strongly favor the presence of NSI 
versus the pure standard interaction. Once the NSI are 
included in the fit, the tension disappears. We find that 
the allowed region in the space of the epsilons (marginal- 
ized over the vacuum parameters) includes the case of 
standard interactions, eee=^eT^^TT— 0, only with confi- 
dence level higher than 99%. 



Note that in both cases shown in Fig. Am^ would 
be measured with approximately 6% precision, and sin^ 6 
with 15% precision. Notice also that, for case (ii), the 
MINOS-high allowed region does not include the input 
values used for the vacuum parameters. This indicates 
that the small corrections due to matter effects on the 
MINOS beam are sufficiently important for the precision 
of MINOS-high. 

We obtain very similar results with a more modest in- 
crease of the MINOS statistics, to 16 • 10^" instead of 
25 • 10^° protons on target. With this intermediate in- 
crease, the point eee=£TT=^eT= is inside the 99% C.L. 
region, but outside the 95% C.L. contour. 
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FIG. 3: Results of fits to simulated MINOS data with statis- 
tics increased from the current 0.93 ■ lO'^" to 25 • 10^" pro- 
tons on target (thin contours). 90% and 99% C.L. regions 
are shown. The "data" were simulate for two sets of NSI 



(i) no NSI, sin^ 
= 0, err= 0.81, e 



0.5 
0.9, 



and "true" oscillation parameters 
and Am^ = 2.7 40"^ eV^ (u) te, 

sin 26* = 0.27 and Am' = 3.1 -10"^ eV^ The fits were done 
in both cases in the assumption of no NSI. For reference, we 
also show the regions allowed currently by all the data com- 
bined, at 90% and 99% C.L. with (filled area) and without 
NSI (thick contours), as in Fig. Q See text for details. 



IV. NSI-DRIVEN i^u 



CONVERSION AT 



MINOS 

So far, we have discussed the dominant — > Vt os- 
cillation mode and argued that the effects of intervening 
matter on this mode are subdominant. We next show 
that for the z^^ — > i/g mode matter effects due to NSI 
can be instead very important. This opens up another, 
complementary mode for searching for NSI at MINOS. 

The z/p Vf, mode is being studied by MINOS as a way 
to measure/constrain the ^13 mixing angle. According to 
the analysis by the collaboration, with the planned higher 
statistics MINOS will have an impressive sensitivity to 
(?i3 in this mode, providing a significant improvement 
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over the current CHOOZ bound. For instance, with 25 x 
2q20 pj-otons on target and assuming Ato23 = 2.5 x lO^'^ 
eV^ the MINOS bound is projected to be sin^ 26*13 < 
0.07, compared to sin^ 20i3 < 0.14 as currently given by 
CHOOZ. 

If the non-standard flavor changing interaction eer is 



present, it will also drive Vn 



conversion. Schemati- 



cally, this conversion can be viewed in two steps: 



l^n > Vt 



(10) 



The first step has already been observed by MINOS, 
with the largest conversion happening in the lower en- 
ergy part of its spectrum (1.5 — 2 GeV). Correspondingly, 
the t^e production according to Eq. (|10|l is also expected 
to peak at low energy. Below, we quantify this by giving 
an approximate analytical expression for the probability 
P(Vji — > and presenting a sample numerical calcula- 
tion. 



V2 neglecting the effects of idttp2 — {^2 — <^i)'02, and 
finally (iii) solve the equation for Vi treating 1/2 obtained 
before as a source, idtipi — <^i2'02- The second step is 
the one containing an approximation, it assumes that 

\{S2-Si)ij2\ > |5l2V'l|- 

It is useful to compare Eq. Hll|l with the exact solution. 



P(2 



y/iS2-Siy+U 



(12) 



This explicitly shows that (|11|) is a good approximation 
so long as ((52 — ^i)^ > 4(5^2- 

Next, consider a three-level problem with 

(0 (5i2 ^13^ 
612 A2 
^13 Ai, 



A. Analytical treatment of i^^ 



conversion 



We begin by noting that the probability Pii^fj, v^) 
will be small (^ 1). This will allow us to treat the process 
perturbatively. To introduce the idea, let us first consider 
a toy example. Suppose we have a two-level problem with 
the Hamiltonian 



H2v — 



Si Si 

S12 (5; 



and the initial state is (0, 1)"^. Assume further that 
1^12! < 1^2 - Si\. Then the probability P(2 1) can 
be computed in perturbation theory as 



P(2^1) 



dt' {-i)5i2e-'^^'-^^^*' 



-i{S2-Si)t 



{S2 - Si) 



(11) 



To arrive at this one can (i) subtract Si on the diagonal 
(overall constants on the diagonal do not change oscil- 
lation probabilities), (ii) solve the evolution equation for 



Notice that the 23 block has a diagonal form. Suppose 
the initial state is i^23 = (0, a, 6)"^ and we again want the 
probability of finding the particle in the first state. If the 
off-diagonal entries in this Hamiltonian are smaller than 
the diagonal ones, we can generalize the previous method 
to give 



P(l^23 



i-i) / dt'iaSue-'^''' +bSi3e-''^'''] 



aSi2- 



bS 



1 - e-^'^i* 



13- 



Ai 



.(13) 



The problem we are interested in, that of Vn 



Ve conversion, can be reduced to the last case. For that, 
we need to diagonalize the — Vr block of the Hamilto- 
nian. This is done by rotating this block by an angle 6*33 
given by 



tan26'n 



A sin 261' 



23 



23 



A cos 26*23 + \/2G'FiVeerr/2' 



(14) 



Notice that this angle is modified from its standard value 
^'23 by the presence of the NSI matter term. 
The two relevant mass splittings are 



Ai,2 ^ i (V2GFNeerr ± ^ {V2GFNeerry + 4(A2 + V2GFNeerrA COS 2^23) ) + A - V2GFNeil + eee), (15) 



sign. Normal hierarchy 
The small solar splitting is 



where Ai comes with the 
is assumed for definiteness 
neglected. 
The state 

sis states; correspondingly, the initial wavefunctions are 



Vfj, is a linear combination of the new ba- 



given by a = cos ^23 and b = — sin ^23 in the rotated 
basis. We also need to rotate the NSI eer term and the 
vacuum 9i2 and ^13 terms. It is important to stress that 
the term is rotated from the flavor basis (by the angle 
6*23), while the other two are rotated from the standard 
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(no NSI) basis that diagonalizes the 23 block (by the an- 
gle - 023). 



J 



Combining these ingredients, we finally obtain 



. , exp(iAiL) - 1 „ „, exp(iA2L) - 1 
Gi sm 6*23 T G2 cos 6I23 V 



Ai 



where 



Gi = \/2Gj^^eeer cos 6*23 + A sin 26*13 cos(6'23-6l23) + A0 sin 26*12 sin(6'23-6l23), 
G2 = \/2Gi^A^eeer sin 6*23 + A sin 26*13 sin(6'23- 6*23) -Aq sin 26*12 cos(6'23- 6*23). 

I 



(16) 



(17) 
(18) 



B. Sample numerical calculation of v^^ 
conversion 




4 6 8 10 12 14 16 18 20 

E,, GeV 



FIG. 4: Conversion probability P(yp. Ve) as a function 
of energy for (i) sin^ 26*13 = 0.07, Amas = 2.5 x 10"^ eV^ 



sin^6*23 = 1/2 
dashed curve), vs. 



and standard neutrino interactions (short- 
) sin^ 26*13 = 0, Am23 = 2.9 x 10"^ 
eV^ sin^ 6*23 = 0.36 and eee = 0, = 0.9, e^^ = 0.81 (solid 
curve). Both curves were computed numerically. For compar- 
ison, the predictions of the analytical result, Eqs. 1161171181 . 
for the case (ii) are also shown (long-dashed curve). 



We are now ready to explore the sensitivity of the — s- 
conversion mode to NSI quantitatively. As already 
mentioned, the MINOS experiment expects to extend the 
lower bound on 6*13 down to sin^ 26* < 0.07. We plot 
the corresponding probability P{v^ Ve) for this case 
in Fig. ^ (short-dashed curve). We also plot the case 
with 6*13 = and values of NSI designed to imitate the 
same effect: fee = 0, — 0.9, err = 0.81 (solid curve). 
Notice that the epsilons were chosen to lie in the presently 
allowed region, by satisfying Eq. ||SJ). The details of the 
parameters are given in the Figure caption. 

The Figure shows that the two cases, while quite dif- 
ferent for <\ GeV, give indistinguishable signals for 
> 1 GeV, the relevant energy range for MINOS. This 
fact, already described in in the context of neutrino 
factories, means that the 3cr bound expected for ^13 could 



be translated into a corresponding bound on the NSI. 

Our analytical results agree very well with those of the 
exact numerical calculations. This can be seen in Fig. 01 
where the long-dashed curve shows Pi^^ Ve) com- 
puted according to Eqs. (|16I17I18|I for the same choice 
of parameters as the solid curve discussed above. The 
agreement between the two is more than adequate. 



C. Discussion: degeneracy of 9ri and NSI effects; 
implications for bounds 

So far we have shown that the NSI are constrained 
under the assumption ^13 = 0. The situation is more 
complicated if the effects of NSI and 6*13 are comparable. 
It is easy to see from Eqs. (I16I17I18|I that the two can 
interfere. The interference can be either constructive or 
destructive, depending on the relative phases of ter and 
e"^ sin 26*. Thus, a non-observation of the con- 
version mode at MINOS strictly speaking would not give 
a simple bound on ^13 or eer, but would define an allowed 
region in the e"' sin 2^13 - e^r parameter space. The de- 
generacy would have to broken by some other means. 



D. Note on MINOS baseline advantage over K2K 

Finally, we note that MINOS, with its three times 
longer baseline, possesses an important advantage in sen- 
sitivity to the NSI relative to K2K. Indeed, consider 
Eqs. (|16I17I18|I and suppose that both the baseline L and 
the neutrino energy Ei, (A oc i/E^) are varied in such 
a way that for shorter baseline the oscillation phase in 
the Vi, v-r channel is preserved. It is then clearly seen 
that the matter term, oc ^/2GFNe€eT, becomes relatively 
less important than the vacuum oscillation term driven 
by 6*13, c>c A sin 2^13. Thus, the K2K experiment is less 
sensitive to the NSI than MINOS. The published K2K 
bound sin^ 26*13 < 0.26 at 90% C.L. |23 does not trans- 
late in any useful bound on the NSI parameter space, as 
we have checked both analytically and numerically. 
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V. SUMMARY 

We have presented the results of a combined anal- 
ysis of the data from atmospheric neutrinos at Su- 
perKamiokande, K2K and MINOS, performed to test 
neutrino nonstandard interactions. The focus has been 
on the effect of adding the recent MINOS data to the 
analysis of the atmospheric and K2K data. 

We find that the allowed region in the space of the pa- 
rameters sin^ 9, Am^, Cee, Crr, Cer h^s the shape predicted 
by the analytics (and confirmed by the data analysis be- 
fore the MINOS data ^3)' has a parabolic direction in 
the space of the NSI couplings and extends to smaller 9 
and larger Am^ with respect to the region found without 
NSI. With the addition of MINOS, the region becomes 
narrower but the part of the region at small 9 due to NSI 
remains allowed, see Fig. ^ 

Another effect of MINOS is to shift the best fit point to 
non-zero NSI, mixing smaller than maximal, and slightly 
larger Arn^. However, there is only an insignificant dif- 
ference in between the best fit point and the point 
that minimizes the for zero NSI. 

Much stronger conclusions could be obtained in the fu- 
ture, when MINOS has higher statistics. Here we have 
shown examples of fits to simulated MINOS data with 
~ 27 times larger statistics ( "MINOS-high" sample) com- 
bined with the atmospheric and K2K data as they are 
today (no updates are expected from K2K, which is now 
closed). We find that in this scenario the precision is suf- 
ficient to reveal oscillation parameters in tension with the 
analysis of atmospheric neutrinos in absence of NSI, and 
thus to give an indication of the existence of new inter- 
actions. It is possible that the point eee=SeT—STT= will 
be outside the 99% C.L. region of the combined analysis 
of all the data. This requires a significant shift of the 
vacuum parameters with respect to the matter ones, 9m 
and Arrim, corresponding to a large "matter" angle, /3, 
/3 ^ O.Stt (see Eq. (|S|). Such condition typically corre- 
sponds to I Cer I ~ 0.5 per electron, i.e. NSI per quark at 
the level of ~ 10% of the Standard Model interaction. 

It should be considered that, since the sensitivity to 
NSI arises from combining data from a neutrino beam 
like MINOS and the signal from atmospheric neutrinos, 
the latter may become the limiting factor as the precision 
of the beam measurements increases. This should moti- 
vate a second phase in the study of atmospheric neutrino 
oscillations, characterized by higher precision, desirably 
of the same level - ~ 10 — 20% - of the perspective MI- 
NOS precision. Such second phase could be realized at 
Megaton water Cherenkov detectors like UNO __28, J9| , 
HyperKamiokande 30] and MEMPHYS Alterna- 
tively, the presence of NSI could be revealed by the com- 
bination of high precision measurements from different 
beam experiments. A typical scenario would require two 
beams of different baseline, one of similar length as MI- 
NOS and the other with a baseline longer by a factor of 
several. For this second beam the matter width would 
be sufficient to have significant refraction effects, if the 



epsilons are close to 1 (per electron, meaning about 
0.2 per u and d quarks). This long baseline experiment 
would play the same role played by atmospheric neutri- 
nos in this paper. Finally, a third possibility to test for 
NSI is to use a MINOS-like beam to look at oscillation 
channels where vacuum oscillations would be suppressed 
by the smallness of the mixing ^13 rather than by the 
high energy, so that small effects of NSI would not be 
hidden by larger vacuum oscillation effects as it happens 
for part of the atmospheric neutrino spectrum. 

We have investigated the latter possibility for the case 
of MINOS-high. The main finding is that epsilons of 
order unity could give a measurable — Ve conversion 
probability of the order of few 10^^. In the interval of 
energy relevant for MINOS the NSI-induced signal for 
^'i3 = can mimic the effect of ^13 (with no NSI), and 
vice versa. To break the degeneracy between the two 
scenarios would require a combination of tests in several 
oscillation channels at a neutrino beam like MINOS itself, 
or a quasi- vacuum test of v^^ — Ve conversion, of the type 
planned for nuclear reactor experiments, e.g. |32|. 

Finally, we again stress the importance of measuring 
the neutral current event rate in the MINOS far detector, 
as a complementary mode for searching for the NSI. Since 
the flux at the far detector contains a significant Vr com- 
ponent, such a measurement would be a valuable probe 
of new tau neutrino interactions. This probe would be 
complementary to the matter effects considered here, be- 
cause the cross-section effects are sensitive to, in general, 
different combinations of the NSI as discussed in (1^] 36l| . 
Such an analysis would be a generalization of the search 
for the sterile neutrino component already planned by 
the MINOS collaboration. 

In conclusion, we have shown that the potential of 
MINOS goes beyond confirming and adding precision to 
the already accepted v^^Vt vacuum oscillations scenario 
that explains atmospheric neutrinos data, and testing 
vacuum oscillations in the Vfj, — Ve channel. This experi- 
ment can be a valuable tool to search for exotic neutrino- 
matter interactions, both in the i'^ disappearance and in 
the Ue appearance parts of its program, and this should 
motivate its long term operation into a phase of pre- 
cision measurements. Our study also makes a general 
point, that it is important to test neutrino oscillations 
with high precision both in environments where matter 
effects are suppressed and in cases with strong matter ef- 
fects, because the comparison of the two could reveal new 
physics. This motivates higher precision measurements 
with atmospheric neutrinos and contributes to make the 
case for long baseline neutrino beams and for a new gen- 
eration of reactor neutrino experiments. 
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onhne, discussing ideas similar to ours about testing NSI 
in the i^^^ — channel at MINOS. Our work adds to 
theirs the three neutrinos analytical formalism, helpful 
to analyze the ^13— degeneracy. 
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